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Ihe Cause of Enches D%msion tn Saturn s Ring. 
By G, E. GoLDSBROUGH, D.Sc, Armstrong College, Newcastle-on-Tjne. 

(Communicated by Prof. T. H. Havelock, F.E.S. Received March 2, 1922.) 

§ 1. Introdttdory . 

In a former paper,* the effect of satellites in producing divisions in Saturn's 
Eings was discussed. The case taken was that where the satellite orbit and 
the rings were co-planar. The results afforded an explanation of the outer 
dimension of the Eing, Cassini's Division, the inner radius of the bright ring, 
and the existence of the Cr^pe Eing. But no reason was given for the 
existence of Encke's Division, nor the numerous divisions reported by Lowell. 
In carrying the examination further, the fact that the satellite orbits are not 
precisely co-planar with the ring system must be considered. As given by 
the ' Annuaire du Bureau des Longitudes,' the inclination of the plane of the 
rings to the ecliptic is 28° 5*6', while the corresponding inclinations of the 
satellite orbits are — 

Mimas 27° 29'6', Dione , 28° 4*4', 

Enceladus 28° 4'3', Ehea 28° 22'8', 

Tethys 28° 40'5^ Titan 27° 39r. 

In the cases of Mimas, Tethys, and Titan, the differences are distinctly 
marked, the first being 36', or O'Ol of a radian. The effect of this inclination 
is examined in the sequel. 

The general hypotheses are the same as in the former paper, to which the 
reader is referred. 

§ 2. The Equations of Motion. 

Using the same notation as in the former paper, let M be the mass of the 
primary supposed situated at the origin of co-ordinates. It is surrounded by 
a single ring of ;p particles, whose unperturbed positions are in the xy plane, 
at the corners of a regular polygon. We shall take the co-ordinates of one of 
these particles, m^, as -X}^, yx> ^a> or r^, d\> Z\ ; r and 6 being measured in the 
xy plane. The unperturbed values of these co-ordinates will be ^a==0, o\^a, 
^^ = ft)^ -f 6 -f- X27r/j), where co^a^ = M -f ?% = M. 

The perturbing satellite, of mass m\ describes a circular orbit of radius a' 
in a plane inclined at an angle 7 to the xy plane. Any ellipticity of the 
primary M is neglected, so that the orbit of m^ has stationary nodes. It 

* * Phil. Trans.,' A, vol. 222, p. 101. 
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has been shown by Stroobant,* that Saturn's Eing, as a whole, has little 
effect upon the movements of the satellites. We are therefore justified in 
treating the orbit of m' as an unperturbed circle. The co-ordinates of m' at 
time i will be 

x' = a' cos d\ y' =z a sin Q' cos 7, ?J = a' sin Q' sin 7, 
where Q' = w'^-f-e', 

and 6)'V^ = M + m^ = M, 

We have here taken the line of nodes as the a?-axis, and have measured Q' in 
the plane of the orbit. 

The equations of motion are now, for the particle mx 



rxQ^ = ; 
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where 
and 



Ax^ = (%-^0' + (^x-/f + (%-^')2, 
= ^\^ + ^/— 2rx7vcos(^x- ^^)+K— %)^. 



In the case of Saturn's Eings, the particles forming any ring can only 
perform small vibrations in the plane of the ring about their zero positions if 
collision with other particles is to be avoided. We therefore assume 

^x = 6^(1 4- px), 
and ^x = coi^ 4- e + \ 2it\]) + a-x, 

where px and <rx are small. 

When the orbit of the perturbing satellite is in the xy plane, %=0 is a 
solution of the equations (1). Further, as we shall only consider the case 
where the inclination of the satellite orbit is small (and therefore we may 
replace sin 7 by 7 and cos 7 by unity), we assume 

fx being also small. 

■^ *Coinptes Rendus/ vol. Vl% p. 913. 
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The equations (1) now become 



a^.+ 2««p. = ^^g^ j^ + tap, (g;~^7^j, 

-. / a^F \ / a^F 

Now 

Ax^ = rx^ + ^A^ + a"^— 2r^a' {cos ^x cos 0' + sin ^x sin 0' cos 7} — 2^x^'. 

Since cos 7 is replaced by unity, and in the formation of zero values of the 
derivatives we put ^^x = 0^ it is clear that, wherever it appears in the 
equations (2), Ax has the same value as in the equations of the former paper. 

In the same wav the other zero values are 

«/ 

Bx^^ = 4a^ sin^ (X — /i) tt/^j, 
x^x' + t/x^/ + ^h^' = aa^ cos (m't 4. c' ^ co^ — € — X ^^Ip)^ 
and a'x^v + yxVt. + z^z^ = a^ cos (X — yu,) 2 Trjjp. 

It is evident, then, that the terms of equations (2) will all be the same as 
those of the former paper except where ^-derivatives apjaear. The latter 
terms are now evaluated :— 

'^'^'^^'^ (a78^)o "^ -^^^|^7'{^^-~^'GOs(e^-~^x)}sin(9'?x 

= 0, to the order of value we are using. 

J'^^" \d^.W.h = A? ''^ ^^ ~^^> ^^ 

= 0, 

''^^■^ = ,;,' sin ^' . Y i -^ - 4- + ^ COS (^' - 0.)| > 



^dzja ' LAa^ a''^ a' 



_, / a^F \ „ , ■ /,/ r aa' — a'^ cos (0' — 01,) ct, ,„' a\\ 



^ X^-^A*^' ^, 



/ 3^r \ r a'^ 1 1 

^^^ ( 5F3~" = ^'^^^' ^''y^^ sin ^' J T-^'sin (^' — ^x) 4- -73 si^^ (^' — ^a) K 



,x W'^M^'^A/0 



A /A -*-'X/x 
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In the last expression the yLt-summation in the right-hand member does 
not include the case />t = X. Also 0' now equals co't + e' and 0x equals 
coif + e + X27r//A 

It is clear from these results that the first pair of equations (2) remain 
unchanged by the small inclination of the satellite orbit. The effect of this 
inclination is only to produce vibrations normal to the plane of the ring. 

We are therefore left with the single set of equations 

C, = m'sin 6' i -^3 h^ + 4s cos(^'~^,) 

4-3m'sm^ ^ ---^ ^4-— cos(0'—^a) [-/^a 

-f3m' sin & ||lsin(^'-^;,) + ~3 sin(r-^x)T <tx 

It is to be noted that this equation holds for all integral values of \ from 
1 to p. It is independent of 7 ; and the quantities pxy cr^ are known for all 
values of \, having been determined from the first two of equations (1) which 
have been solved to the necessary degree of approximation previously. 

In the former paper, certain approximations were made to render possible 
the solution of the period equation without a definite knowledge of the 
separate masses forming a ring.* This approximation amounted to the use of 
the mean position instead of the actual position of a particle in calculating 
its perturbing effect. upon another particle. If we adopt the same process 

again, we reject the terms 2 fr^ Tm i^ equation (3). We have then a series 

of p equations for the motion of the p particles, each equation containing the 

co-ordinates of the particle under consideration and not the co-ordinates of 

the other particles. This of course very much simpHfies the work. 

We may write 

Ax"2 = a'-2Sc^cos^((9'-6>A), 

i 

Aa"^ = a' -'' 1, di COS i{0'~ did. 

i 

On substituting these and reducing the expressions, we find that equations (3) 
may be written 

r, + ?,2<i>i,,-cosy(0'-^;,)-fp,S<i>2,,-sin{o + i)0'-y^,} 

-^a,X^,,jOOB{(j + l)0'^je,} = 2^4„sin {(j+l)0' ^j 0,}, (4> 
^ Log. cit.^ pp. 123 and 124. 
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In these equations the vakies of the coefficients $ are readily obtained. In 
particular 



hd 



m"" -tco'^ 2^ -YF -T-Yjz- r— T~ -f-V CD -Co. 



Ab was shown in the previous paper * % -z— -~-t--— ~ — — rans^es in 

^ ^ ^ ^ M 8 sm^ (X — m) tt/p ^ 

value from zei^o to 0*0096 y^^m/M, where m is the maximum mass of a particle 
appearing in the ring. All the quantities # except #1,0 are factored by v' or 
m'/M. As this quantity is very small, all the quantities # except ^1,0 will 
be small, and their squares and products may be neglected. Also in #1,0 the 
term v'q}'% may be rejected in comparison with the other two. We shall 
then write 

§ 3. The Solution of the Equations, 

In solving equations (4) let us first take the complementary function. The 
equation to be solved is 

?+ f [2#i,iCos/(€a'^ + e^-ft>^-e--.X-2'?r/p)] = 0. (o) 

The solutions of the system are obtained by giving to \ its j? values 

1, 2, ...^. The equation is of HilFs well-known form and can readily be 

solved by the method of Whittaker and Ince.f To simplify the form, change 

the independent variable from t to ^ which equals w'^ + e'— ct)^— e— X27r/p. 

Then 

1 



r+r 



]S#i,jCosy^ =0. (6) 



Let ^ = e^hi, where u is a periodic function of the independent variable ^. 



Then 

1 



Also let 



S#i,jcos/^ 



0. (7) 



i 

f ~ VO<I>, • 

i 

In the summations, all integral values of i are to be taken except ^ = 0, 
The general case should include products of the quantities #, but as has been 
pointed out these may be neglected. The quantities Q^ and E*- are constants, 

* P. 123, 

t « Monthly Notices, R. A.S./ LXXV, 5, p. 436. 
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and Pf is a periodic function of <^. If we add the two conditions : (1) That 
no term mx(n<^—r) may appear in the expression for ze, and (2) that 'i^ is 
purely periodic, we are able to determi^ie completely all the coefficients 
entering into the expressions. 

On substituting in equations (7) and equating to zero the coefficients of 
those terms which are independent of #, we, have an identity. From the 
coefficient of #i, i we find 

- 2n Q,- sin {nd> - r) + V-' + n^l^i + , ^' ^. cos (n6 - r) 

(ft) — COf 

, cos(t^c(>— T)cos'i<i A /Qv 

In general, therefore, 

Q^ = 0, E,- = 0, 

p, _ -hG0s{(i-hn)4>—T} i cos {(i—n) if) --t} xq>. 

* (ft>'-ft>)2 {(iJ^-nf-^n^yiay'-^-wf {(i^-nf-n^Y ^ ^ 

But in the particular case of i = 2% equation (8) reduces to 
— 2n Q2n sin (ncj} - t) + Psn' ' + n^ V2n + , ^^"^ , , cos (?i6 - r) 

(ft) —coy 
1 

+ ~--j—-~ {I cos(3^<)b — T) 4- 1 cos(^i^--t)cos 2t--| sin (?i<^ — r) sin 2t} = 0. 

Whence, to fulfil the given conditions, 

sin2T 



We have, then, 



4:71 {(O -^wf 

E2» =: -~|cos2t, 

p _ cos<3^— t) ,.^, 

'"^ '^ 16^^(ft)'-ft)f ^^^^ 

472,^(0 — ft))^ 

^and #1,0 = #(c»'-ft>)2-|#i.2„cos 2t. (11) 

These two expressions, together with (9) and (10), constitute a solution 
valid in the vicinity of *i,o = ^i2(ft)'~ft))2, where n is any integer. The 
parameter t is determined from the second of equations (11). The value so 
obtained, substituted in the expression for c, determines that quantity 
completely. It is clear that real values of r mean real values of c, and hence 
vibrations of increasing amplitude or instabihty. Real values of r can only 
occur when 

|<E>i,o-^i^(ft)'-ft))2[^||<l>i,2.ii ; 
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and owing to the extreme smallness of <l>i, 2m this is, with a close degree of 

approximation, when 

*i,o-~^'K-a))^^0. (12) 

We pass next to the particular integral of equation (4). The values of p}^ 
and <7x to be used are obtained from the former paper. From p. 110, we 

have 

p = e^^ {Ao sin (n^—r) + %lAr, s @r, « . • . }, 

and a = ^^^ {Xocos(w<^— T) + S2Xy,5%,s -. }• 

In these expressions Ao and Xo are related constants^ and A^, s, X^, s are 

determined later in the paper. As, however 0^,^ is of the same order as 

€>2, j and <I>3,j, it is evident that on substituting in equation (4) only the first 

terms of the expressions for p and a- need be used. We have further to add 

the particular integrals for p and a given on p. 117 of the former paper. As 

again, these particular integrals are of the same order of value as the 

quantities ^2,j and #3,j it is clear that we may omit them in making the 

substitutions. The final result is that the values of p and cr will both produce 

terms of the type 

Cli COB {(i+l) 6' — id],} . GX^.c{0' — Ox), 

where i may take all integral values, positive and negative, and c and O^ are 
known constants of the same order as ^2>j and ^s,j. 

Particular integrals will arise from these terms and also from the remaining 
terms 

^,,jsin{{j+i)d'--je,}. 

Taking the first of these we have to find the particular integral of 

r,+ rAS^i>iCOs/(r-^x) = a-cos {{i+lje'-^idx} exp.c(r-0x). (13) 

j 

Eemoving exp. {e (^'•— ^a)}> and neglecting c^ we have 

J 
On working oat the particular integral in the way previously adopted,* and 
noting that products of # and O may be rejected, we have 

or, to the same order of values, rejecting the product cili, 

r,= ^"^P'i^j^\-"/^)^%cos{(^ + l)g--^gA}. (14) 

^ Loc. cit., p. 116. 
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To this must be added the particular integral of 

To the same order of approximation, the required integral is 

>- _ <a>4,,-siii{o-+i)g^-ygx} 

-^ *i,o-{0'+l)«'-ya)}2' 



(15) 



(16) 



Finally, the whole particular integral is represented by the sum of (14) 
and (16). 

§ 4. Discussion of tfie Results, 

Expression (12), with the meaning of ^i,o inserted, is 

a)2-fft>2F;, = #(c«'~a>)2. (17) 

As already pointed out, for the differing values of X, F^ ranges from zero to 
0*0096 j9^m/M. As in the previous paper, we shall treat F^ as a parameter 
varying continuously between the limits mentioned. Solving (17) and taking 
only the case where &>>«', 

^' = ^-\/(l + F;,) ^ 

(o n ' ' 

The maximum value of F^ as determined from the observed dimensions of 
Cassinfs Division is 0*0173. The following Table gives the positions of the 
divisions resulting from the instability : — 





JFa - 0. 


Pa = -oils. 


od'/oo. 


a/a\ 


(v'/co. 


aja'. 


2 
3 



i 

3 




0-630 

0-763 


-4958 
-6639 


'626 ■ 
0-761 



Applying these results to the case of satellite Mimas, we should find, from 
the solution ^2. = 2, a narrow division in Saturn's Eing extending from 16-9" 
to 16*8". This is just at the inner edge of Cassini's Division, For Enceladus 
the same solution would give a narrow division from 21*7'' to 21-6''. This is 
outside the limits of the existing ring. The solution for i^i = 3 gives for 
Mimas a narrow division at 20'5". This is again just outside the existing ring. 

The normal vibrations represented by the solution (17) do not lead, then, to 
any new divisions in the Eing of Saturn. The particular integral, represented 
by the sum of (14) and (16), becomes large, and instability ensues in the 
vicinity of that position where 

<3E>i,o-{(^ + l)ft)'-iei>}2 ~ 0. (19) 
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Taking the value of ^i^o, equation (19) becomes 

ft)2(l + Fx)-{(^+l)6>'-~^(oP = 
The solution of this equation for which &) > «' is 



0. 



ft) 
ft) 



i + 1 



(20) 



As before, Fx ranges in value from to 0*0173. The following Table gives 
the limits of the zones of instability arising from (20) for various values of i, 



which is an integer :- 





FA = 0. 


1 
Fa - -0173. 


1 

co7(w. ■ 


a/a\ 


(a'jo). 


aja' . 


i = 1 

2 
3 

4 
5 




1 

3 
JL 

a 

3 
5 
3 

3 




0-481 

0-630 

0-711 

0-763 


0-331 
-498 
0-598 
0-665 


0-478 
-628 
0-710 
0-762 



As appears in this Table, a series of narrow divisions arise from this part 
of the integral. 

On applying the ratios to the satellite Mimas, we iiave the following series 
of divisions. They are worked to two places of decimals, in order to show the 
order of magnitude of the width of the divisions. 

For n = 2, 12-90'' to 12*81'^ 
^ := 3, 16-89'' to 16-85", 
^ = 4, 19-08" to 19-03", 
n == 5, 20-50" to 20*45". 

The first division is in. the Crepe Eing, just within the commencement of 
Eing B, the inner diameter of which is 13-00". 

The second division is at the inner edge of Cassini's Division. 

The third division is a new one and falls almost exactly in the position 
given for Encke's Division. Lowell* gives the radius of Encke's Division as 
19-00" ; while Seef states it at 18*87" with a width of 0*41". In either case 
theory clearly indicates the existence of this division. 

The third division again falls outside the existing ring. 

Applying next the ratios in the Table to the satellite Enceladus, the only 

case to be mentioned is the value when % = 2. This should give a division 

extending from 16-56" to 16-44". 

^ * Lowell Observatory Bulletin,' No. 68. 
t ' Astr. Nach.,' No. 3768. 
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There is no observed division at this position that I have been able to find. 
But it is easy to show that it must be very indefinite. The value of the 
inclination of the orbit of Enceladus is 1'3' compared with 36' for that of 
Mimas. The effect of the inclination must, therefore, be less marked in the 
former case. Again, if we compare the distance of the satellite producing 
the Encke Division from that division with that of Enceladus from the new 
division, taking ratios (as that is the form in which the distances appear in 
the analysis) we find a/a^ = 0**7 for Encke's Division and a/a^ = 0*48 for the 
new division. The perturbing effect of Mimas in Encke's Division is, therefore,, 
much greater than that of Enceladus in the new division. On these two 
counts we should expect the new division to be less distinct than Encke's 
Division. But Encke's Division is only visible with difficulty ; hence it is to 
be expected that the new division should be unobservable. 

None of the remaining satellites produce divisions corresponding to the 
Ibrmulse of this paper within the existing Eing system. 

I 5. Summary, 

The effect of a satellite moving in an orbit inclined to the plane of the 
Eings of Saturn is to produce one new division. If the satellite be taken as 
Mimas, a narrow division will be caused at a distance 19*05'' from Saturn. 
This closely corresponds with the observations of Encke's Division. 

Enceladus should produce a division at 16'5". No observed division is 
recorded at this place, but because of the minute inclination of the orbit of 
Enceladus and because of this satellite's distance from this place, it is to be 
expected that the division will be comparatively feeble. 



